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ABSTRACT 

The  time-averaged  intensity  density  function  of  the  acoustic 
radiation  from  rotating  blades  is  derived  by  replacing  blades  with 
rotating  dipoles.  This  derivation  is  done  under  the  following  turbu¬ 
lent  inflow  conditions:  turbulent  ingestion  with  no  inlet  strut 
wakes,  inflow  turbulence  elongation  and  contraction  with  no  inlet 
strut  wakes,  and  inlet  strut  wakes. 

Dimensional  analysis  reveals  two  non-dimensional  parameters 
which  play  important  roles  in  generating  the  blade-passing  frequency 
tone  and  its  multiples.  The  elongation  and  contraction  of  inflow 
turbulence  has  a  strong  effect  on  the  generation  of  the  blade-passing 
frequency  tone  and  its  multiples.  Increasing  the  number  of  rotor 
blades  widens  the  peak  at  the  blade-passing  frequency  and  its  multi¬ 
ples.  Increasing  the  rotational  speed  widens  the  peak  under  the 
condition  Chat  the  non-dimensional  parameter  involving  the  rotational 
speed  is  fixed.  The  number  of  struts  and  blades  should  be  chosen  so 
that  (the  least  common  multiple  of  them) • (rotational  speed)  is  in  the 
cutoff  range  of  Sears'  function,  in  order  to  minimize  the  effect  of 
the  mean  flow  deficit  on  the  time  averaged  intensity  density  function. 
The  acoustic  intensity  density  function  is  not  necessarily  stationary 
even  if  the  inflow  turbulence  is  homogeneous  and  isotropic.  The  time 
variation  of  the  propagation  path  due  to  the  rotation  should  be- con¬ 
sidered  in  the  computation  of  the  intensity  density  function; 
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for  instance,  in  the  present  rotor  specification,  the  rotor  radius  is 
about  0.3  m  and  the  rotational  speed  Mach  number  is  about  0.2. 
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CHAPTER  I 


INTRODUCTION 

1.1.  Statement  of  the  Problem 

This  scudy  attempts  to  derive  the  time  averaged  far  field 
acoustic  intensity  spectral  density  function  due  to  the  acoustic 
radiation  from  (1)  low  rotational  speed  rotors  operating  in  frozen, 
homogeneous,  and  isotropic  turbulence,  and  (2)  low  rotational  speed 
rotors  behind  inlet  struts,  where  the  rotor  is  assumed  to  be  acousti¬ 
cally  compact.  For  the  present  analysis,  the  dipole  source  due  to 
blade/turbulence  interaction,  as  computed  by  Sears'  function,  is 
assumed  to  be  predominant  over  other  radiation  sources.  A  rotating 
dipole  model  is  pursued  so  that  the  effect  of  the  dipole  rotation  can 
be  accounted  for  explicitly  in  the  analysis.  Sears'  function  is 
regarded  as  an  impulse  response  function  of  a  blade,  so  that  the 
analysis  is  done  in  a  time  domain  and  in  a  frequency  domain.  There 
is  no  blade-to-blade  interaction  to  be  included,  such  as  that  devel¬ 
oped  by  Kemp,  et  al.  [1].  It  is  assumed  that  the  inlet  struts  do  not 
affect  the  sound  propagation  and  that  the  distance  from  the  rotor  to 
the  duct  inlet  or  exit  is  much  shorter  than  the  sound  wavelength  so 
that  duct  effects  on  propagation  are  negligible. 

1.2.  Previous  Investigations  of  Rotor  Noise 

Rigorous  treatments  of  the  effect  of  a  solid  boundarv  under  the 
influence  of  a  volume  distribution  of  quadrupoles,  as  in  the  Lighthill 


acoustic  analogy  approach  [2],  were  done  by  Curie  [3]  and  Doak  [4] 
for  a  stationary  boundary,  and  by  rfowcs  Williams,  et  al .  [5]  for  a 
non-stationarv  boundary  and  a  rigid  surface  boundary  of  a  moving  bodv 
such  as  that  of  a  rotating  blade  with  a  low  aspect  ratio.  Curie 
derived  the  sound  field  as  the  sum  ot  that  generated  by  a  volume 
distribution  of  quadrupoles  and  by  a  surface  distribution  of  dipoles 
due  to  the  diffracted  wave  from  the  quadrupole  and  due  to  the  hydro/ 
aerodynamic  flow  itself;  but  the  distribution  of  the  dipole,  such  as 
the  tangential  component  of  the  dipole  on  the  surface,  is  difficult 
to  measure.  Doak  obtained  the  sound  field  for  the  same  situation  as 
did  Curie,  but  the  former's  result  concerning  the  contribution  from 
the  surface  dipole  is  based  on  an  easily  measurable  quantity.  Ffowcs 
Williams,  et  al.  obtained  the  sound  field  as  the  sum  of  that  generated 
by  (1)  the  quadrupole  in  a  volume  distribution,  (2)  the  dipole  in  a 
surface  distribution,  and  (3)  the  monopole  in  a  surface  distribution 
where  the  monopole  is  due  to  the  volume  displaced  by  a  non-stationary 
boundary.  He  found  also  that  in  case  (3)  as  for  the  case  of  the  rigid 
surface  boundary  of  the  moving  body,  the  monopole  is  broken  into 
another  quadrupole  and  dipole. 

In  the  Lighthill  acoustic  analogy  approach  [2],  the  acoustic 
radiation  source  region  is  assumed  to  be  substantially  small,  because 
the  observer  of  the  radiation  is  far  from  the  source  region,  so  that 
the  propagation  and  source  terms  are  separable  and  the  non-linear 
effect  is  regarded  as  the  source  term  in  the  acoustic  wave  equation. 

In  other  words,  the  source  distribution  does  not  affect  the  acoustic 
propagation  or,  more  preciselv,  the  region  where  the  acoustic  field 
is  to  be  computed  does  not  overlap  with  the  source  region,  so  that  no 
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convection  effects  are  involved  in  the  acoustic  wave  propagation. 

Thus,  if  the  linear  acoustic  process  is  isentropic,  then  the  acoustic 
intensity  is  computed  based  on  the  acoustic  pressure. 

In  the  Lighthill  acoustic  analogy  approach,  the  existence  of  a 
rigid  body  yields  dipole  and  quadrupole  terms  in  the  acoustic  wave 
equation.  However,  their  relative  importance  is  dependent  on  the 
rotor  speed,  the  dimensional  parameters  of  the  rotor,  and  the  fre¬ 
quency  range  of  interest.  In  practice,  there  are  situations  where  the 
dipole  source  dominates  the  acoustic  radiation  over  that  due  to  the 
quadrupole  source.  Ffowcs  Williams,  et  al.  [5]  and  Goldstein  [6] 
showed  the  relative  importance  of  the  dipole  source  over  the  quadru¬ 
pole  source  for  a  low  speed  rotor  with  a  small  number  of  blades,  under 
the  assumption  that  the  dipole  and  quadrupole  sources  have  similar 
spatial  and  temporal  scales. 

The  assumption  of  acoustically  compact  blades  permits  replacing 
distributed  dipole  sources  with  a  single  point  dipole  whose  strength 
equals  the  total  strength  of  the  distributed  dipoles  on  the  blade. 
Thus,  each  blade  is  replaced  with  a  point  dipole  or  multiples  of  a 
dipole.  Lowson  [7]  derived  the  far  field  and  near  field  acoustic 
pressure  from  a  moving  singular  dipole  source  by  using  Lighthill's 
acoustic  analogy  approach,  where  the  importance  of  the  acceleration  of 
the  moving  dipole  was  accounted  for  in  the  expression  of  acoustic 
pressure . 

As  there  are  two  interpretations  of  the  solution  of  the  acoustic 
wave  equation  of  cylindrical  coordinates,  i.e.,  a  rotating  pressure 
pattern  using  an  exponential  function  for  the  angular  variation  of  the 
pressure  pattern,  or  a  modal  pressure  pattern  fluctuating  in  time 
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using  a  trigonometric  function  tor  the  angular  variation  of  the 
pressure  pattern,  the  acoustically  compact  rotating  blades  can  then 
be  modeled  by  on-off  dipoles  or  rotating  dipoles.  Horaicz  [8],  Gutin 
[9],  Wright  [10,  11],  and  others  pursued  the  on-off  dipole  model. 

Ffowcs  Williams,  et  al.  [12],  Lowson,  et  al.  [13],  Morfey,  et  al.  [14], 
Mani  [15],  and  others  employed  the  rotating  dipole  model. 

The  rotating  dipole  model  accounts,  on  the  one  hand,  for  the 
fact  that  the  distance  between  the  observer  and  the  rotating  dipole 
varies  according  to  tne  rotation  of  the  dipole,  and  that  the  direction 
of  the  motion  of  the  dipole  varies.  On  the  other  hand,  the  on-off 
model  does  not  account  for  it  because  the  distance  between  the  on-off 
dipole  and  the  observer  is  fixed  all  the  time. 

There  are  two  ways  to  get  the  dipole  strength  necessary  to 
compute  the  acoustic  radiation  from  rotating  blades:  analytically  and 
empirically.  For  instance,  Hanson  [16]  derived  the  far  field  acoustic 
pressure  for  a  compact  subsonic  rotor  based  on  rotating  dipoles  whose 
strength  was  derived  from  the  pressure  measured  by  a  pressure  sensi¬ 
tive  transducer  mounted  at  a  specific  point  on  the  blade  surface,  even 
chough  the  position  of  the  transducer  has  a  significant  influence  on 
Che  acoustic  pressure  ac  the  observer's  point.  In  addition,  the  size 
of  the  pressure  transducer,  i.e.,  the  area  of  the  transducer,  influ¬ 
ences  the  measured  data  especially  ac  high  frequencies.  Homicz  [8], 
Aniec  [17],  Aravamudan,  et  al.  [IS],  Mani  [19],  Sevik  [20],  and  other 
utilized  a  two-dimensional  aero/hydrodynamic  function  developed  by 

Sears  [21]  in  order  to  obtain  the  dipole  strength. 

% 

The  two-dimensional  aero 'hydrodynamic  function  was  developed  for 


a  two-dimensional  thin  airfoil  having  small  camber  under  the  following 


assumptions:  (1)  the  flow  from  the  upper  and  lower  surfaces  of  the 

blade  at  the  trailing  edge  has  no  pressure  jump,  (2)  a  thin  sheet  of 
vorticity  is  shed  from  the  sharp  trailing  edge  (Kutta  condition), 

(3)  the  thin  sheet  is  shed  along  the  chord  line  of  the  blade,  (1)  the 
finite  thickness  of  the  boundary  layer  at  the  trailing  edge  is  not 
important,  (5)  no  viscosity  is  involved  although  the  vorticity  is 
generated  by  viscous  action  in  the  airfoil  boundary  laver,  (6)  gusts 
to  the  airfoil  are  frozen  (the  frozen  gust  assumption  implies  that 
gusts  are  not  distorted  during  the  time  of  passage  over  the  airfoil) , 
and  (7)  the  airfoil  is  acoustically  compact. 

Among  the  aforementioned  assumptions,  the  Kutta,  frozen  gust, 
and  compactness  conditions  require  special  care  when  using  the  two- 
dimensional  aero/hydrodynamic  function  for  the  following  reasons: 
although  the  Kutta  condition  has  a  dramatic  effect  on  the  pressure 
distribution  on  the  blade  surface,  as  shown  by  Kelly  [22],  the  condi¬ 
tion  is  still  being  investigated  [23].  The  frozen  gust  condition 
implies  that  the  two-dimensional  aero/hydrodvnamic  function  is  rela¬ 
tively  valid  for  low  values  of  the  frequency  compared  with  higher 
values.  The  assumption  about  compactness  asserts  that  the  acoustic 
radiation  at  one  surface  point  does  not  affect  the  surface  pressure 
at  others.  Thus,  at  high  frequencies,  it  is  meaningless  to  integrate 
the  surface  pressure  in  order  to  get  the  strength  of  a  point  dipole. 

Jackson,  et  al.  [24]  compared  the  measured  and  calculated  aero¬ 
dynamic  admittances  which  are  associated  with  aero/hydrodvnamic 
functions  such  as  are  given  by  Sears'  two-dimensional  theorv,  bv  two- 
dimensional  strip  theory,  and  by  three-dimensional  theory.  The  two- 
dimensional  strip  theory  states  that  the  lift  of  each  chordwise  wing 


element  is  a  function  only  of  the  local  upwash  fluctuations  at  the 
element,  and  that  the  lift  on  each  strip,  i.e,,  each  element,  is 
obtained  by  using  Sears'  two-dimensional  theory.  Results  show  that 
the  three-dimensional  theory  such  as  that  of  Filotas  [25]  or  Graham 
[26]  gave  good  agreement  with  measurements.  The  two-dimensional  strip 
theory  gave  gairly  good  agreement.  Sears’  two-dimensional  theory 
overestimated  for  the  entire  reduced  frequency  range. 

Amiet  [27]  showed  that  the  two-dimensional  strip  theory  is  a 

good  approximation  for  computing  the  far  field  acoustical  pressure  due 

to  unsteady  lift  if  >tk  d  is  large,  where  M,  k  ,  and  d  are  a  streamwise 

x  x 

Mach  number,  a  chordwise  turbulence  wave  number,  and  a  span  length, 
respectively.  The  most  important  development  in  his  argument  is  that 
the  computed  lift  agrees  quite  well  with  the  measured  lift  even  when 
Mk^d  =  1.  He  showed  that  the  acoustic  radiation  produced  by  the 
parallel  gusts  would  be  the  most  dominant  noise  source. 

The  two-dimensional  aero/hvdrodynamic  function  causes  aliasing 
among  different  wave  number  gusts.  This  is  because  the  two- 
dimensional  aero/hydrodynamic  function,  such  as  Sears'  function, 
utilized  a  wave  number  in  the  direction  of  the  chord  line  of  a  thin 
airfoil.  Therefore,  even  if  the  frozen  gust  is  decomposed  into  its 
wave  number  components,  each  component  cannot  be  the  same  among 
rotating  blades,  because  the  direction  of  the  chord  line  of  each 
blade  is  different  from  each  other. 

It  should  be  noted  that,  since  the  two-dimensional  aero/ 
hydrodynamic  function,  such  as  Sears'  function,  involves  one  reduced 
frequency  relating  the  wave  number  of  the  upwash  to  an  airfoil  to  a 
frequency,  then  the  interaction  between  blades  and  the  flow 
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fluctuation  can  be  treated  in  a  title  domain  and  in  a  frequency  domain 
instead  of  a  space  domain  and  a  wave  number  domain. 

For  the  analysis  of  the  acoustic  radiation  from  rotating  blades, 
there  are  three  possible  flow  types:  (1)  steady  uniform,  (2)  steady 
non-uniform,  and  (3)  unsteady.  Gutin  [9]  investigated  the  acoustic 
radiation  from  rotating  blades  under  a  steady  uniform  flow,  i.e.,  a 
constant  lift  force.  Morse,  et  al.  [ 23 ] ,  Lowson,  et  al.  [13],  and 
Wright  [10,  11]  tried  to  explain  the  discrepancy  between  Gutin's 
computed  results  and  experiments,  especially  in  the  higher  order 
harmonics  of  the  blade-passing  frequency,  by  introducing  a  steady  non- 
uniform  flow.  As  quoted  by  Lowson  [29],  Leverton  [30]  reported  that 
no  blade-passing  frequency  and  its  harmonics  appeared  in  the  spectrum 
of  the  acoustic  radiation  from  a  rotor  under  calm  wind  conditions. 
However,  a  slight  change  in  conditions  caused  a  blade-passing  fre¬ 
quency  and  its  harmonics.  In  a  steady  uniform  and  a  steady  non- 
uniform  flow,  the  phase  difference  of  the  loading  on  different  blades 
can  be  expressed  deterministically  by  the  geometrical  angle  position 
difference  among  blades.  This  is  not  the  case  for  an  unsteady  flow. 

Homicz  [3]  took  the  on-off  dipole  model  as  the  representation  of 
rotating  blades  operating  in  an  unsteady  flow  of  convected  frozen 
turbulence  by  replacing  a  blade  with  a  dipole  at  a  representative 
radial  position.  Sears'  function  was  used  to  compute  the  lift  on  a 
blade,  where  a  reduced  frequency  averaged  over  one  blade  revolution 
was  used  in  order  to  get  around  the  aforementioned  aliasing  phenom¬ 
enon.  He  showed  that  the  total  power  of  the  acoustic  radiation  from 
rotating  blades  operating  in  homogeneous  and  frozen  turbulence  is 
proportional  to  w~U  where  w  and  l'  are  a  mean  square  turbulence 
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velocity  and  a  convection  velocity  of  the  turbulence,  respectively. 

He  showed  further  that  the  blade-passing  frequency  and  its  multiples 
occur  in  the  frequency  spectrum  of  the  acoustic  radiation  when  U/(Aft) 
is  very  small  where  U,  ft,  and  A  are  an  axial  flow  velocity,  a  rota¬ 
tional  speed,  and  a  typical  length  scale  of  the  turbulence,  respec¬ 
tively.  The  smallness  of  U/(Aft)  implies  chat  a  partially  coherent 
turbulence  is  chopped  by  many  blades  so  that  blade-passing  frequency 
and  its  multiples  appear  in  the  frequency  spectrum  signature. 

Mani  [15]  investigated  the  noise  due  to  the  interaction  of  inlet 
turbulence  with  isolated  stators  and  rotors.  He  assumed  that  the 
amplitude  strength  for  each  rotating  dipole  is  equal  and  did  not 
consider  the  phase  relation  of  the  strength  of  dipoles,  which  should 
occur  for  turbulence  interaction.  The  blade-passing  frequency  and  its 
multiples  in  the  frequency  domain  were  assumed  due  Co  the  aerodynamic 
interference  between  moving  blade  rows  as  investigated  by  Kemp,  et  al. 
[1].  He  showed  that  the  sound  power  level  increases  with  decreasing 
L/D  where  the  turbulence  is  assumed  to  be  homogeneous,  isotropic,  and 
stationary,  and  can  be  characterized  by  a  longitudinal  correlation 
function  such  as  exp(-r/L)  where  L  is  the  length  scale  of  the  turbu¬ 
lence,  and  D  is  a  transverse  spacing  between  blades.  He  also  showed 
that  the  acoustic  spectra  for  rotors  are  peaked  at  the  blade-passing 
frequency  and  its  multiples  when  L/D  exceeds  about  0.5,  and  chat  when 
L/D  <0.5  the  spectra  start  broadening  noticeably.  However,  the 
parameter  L/D  does  not  take  into  consideration  the  contribution  of  the 
rocational  speed  of  the  rotor  and  the  axial  flow  velocity.  That  is, 
even  if  L  were  small,  a  higher  rotational  speed  and  lower  axial  flow 
velocity  would  create  Che  same  situation  as  when  the  length  scale  of 
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turbulence  is  large,  together  with  a  high  axial  flow  velocity  and 
lower  rotational  speed. 

Sevik  [31]  derived  the  sound  power  spectrum  from  a  rotor  operat¬ 
ing  in  turbulence.  The  noise  source  is  the  axial  component  of  the 
blade  lift  only.  The  sound  power  is  computed  based  on  the  general 
relation  S  (w)  =  E{ ] X(w) J where  S  )  is  the  power  spectrum  of  a 

30 

stochastic  process  x(t)  satisfying  j  otR  (a)da  <  where  R  (a)  is 

XX  XX 

.00 

the  correlation  function  of  x(t)  which  is  assumed  to  be  a  wide-sense 
stationary  process,  and  the  correlation  function  does  not  account  for 
the  rotation  of  the  blade.  The  assumption  of  the  correlation  function 
regards  the  rotor  as  being  stationary.  3ecause  of  the  use  of  only  the 
axial  component  of  the  lift,  the  result  applies  only  near  the  center- 
line  axis  of  the  rotor.  In  Sevik' s  calculation  of  the  lift  of  a 
blade,  the  blade  is  replaced  by  a  line  dipole  source  and  the  idea  of 
a  typical  section  is  introduced,  whereby  the  resultant  velocity  and 
chord  of  the  various  elements  (or  infinitesimal  strips)  may  be 
represented  by  those  of  a  single  "typical  section"  located  at  some 
fraction  of  the  span  of  a  blade.  Unfortunately,  the  resulting  power 
spectrum  does  not  predict  the  blade-passing  frequency  component  and 
its  multiples,  due  to  the  assumption  of  the  correlation  function  of 
the  upwash  at  the  blades,  although,  as  shown  by  Robbins,  et  al.  [32], 
the  level  and  general  shape  of  the  predicted  spectra  agree  with  the 
experimental  data. 

1.3.  Scone  of  Investigation 


In  Chapter  II,  an  acoustic  radiation  formula  is  obtained  from 
the  inhomogeneous  acoustic  wave  equation  with  rotating  dipole  sources. 
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This  leads  to  the  integral  equation,  with  Green's  function,  of  the 
acoustic  wave  equation  in  an  infinite  medium.  The  acoustic  compact¬ 
ness  condition  is  introduced,  so  that  the  blade  can  be  replaced  by  a 
point  dipole.  The  time-averaged  acoustic  intensity  spectral  density 
function  is  developed  in  a  general  form  where  the  effect  of  the 
rotation  of  the  dipole  is  shown  explicitly. 

In  Chapter  III,  the  impulse  response  function  of  a  blade  based 
on  two-dimensional,  unsteady,  thin  airfoil  theory  is  discussed  so  that 
the  blade  response  to  turbulence  can  be  treated  in  a  time  and  fre¬ 
quency  domain.  In  addition,  the  correlation  functions  of  che  lift  on 
blades  and  the  lift  derivative  with  respect  to  time  are  obcained. 

In  Chapter  IV,  the  power  spectrum  of  the  upwash  to  the  blades  is 
obtained  in  the  following  cases:  (1)  the  inflow  turbulence  is  homo¬ 
geneous,  isotropic,  and  frozen,  (2)  the  turbulence  in  (1)  is  elongated 
and  contracted  on  its  course  to  the  rotor  plane,  and  (3)  there  are 
wakes  due  to  inlet  struts  in  addition  to  the  inflow  turbulence.  The 
results  of  Lane  [33]  are  utilized  in  the  derivation. 

In  Chapter  V,  the  results  in  Chapter  II,  III,  and  IV  are 
combined  in  order  to  obtain  the  time-averaged  acoustic  intensity 
spectrum.  The  condition  for  neglecting  the  retarded  time  difference 
among  the  lift  of  blades  is  introduced.  Furthermore,  important  non- 
dimensional  parameters  are  introduced  so  that  the  effect  of  the 
various  turbulence  properties  and  che  rotor  geometry  and  operating 
condition  can  be  determined. 

In  Chapcer  VI,  computations  are  done  according  to  the  results  in 
Chapter  V  such  as  Eg.  (5.1b)  and  Eq.  (5.22)  by  changing  the  number  of 


struts  and  rotors,  Che  spacing  between  struts  and  a  rotor,  che 
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observation  angle,  and  non-dinensional  parameters  such  as  and 

R_/.'w.  In  particular,  the  variation  of  these  non-dimensional  param¬ 
eters  are  made  around  an  axial  flow  velocity  of  =  30  m/sec,  a  rotor 
radius  of  R,  =  0.30  m,  a  rotor  rotational  speed  of  2  =  220  rad/sec. 


an  inflow  turbulence  length  scale  of  A^  =  0.15  m  and  0.017  m,  and  a 
wake  turbulence  length  scale  A  =  0.0063  m. 


CHAPTER  II 


ACOUSTIC  RADIATION  FROM  ROTATING  DIPOLE  SOURCES 

2.1.  Introduction 

There  are  two  tasks  necessary  for  analyzing  the  acoustic  radi¬ 
ation  from  rotating  blades:  (1)  obtaining  the  lift  on  the  blades  due 
to  the  interaction  between  the  rotating  blades  and  turbulence  in  a 
flow,  and  (2)  obtaining  the  radiated  sound  from  the  lift.  In  this 
chapter,  the  acoustic  radiation  due  to  lift  on  the  blades  is  consid¬ 
ered,  based  on  the  Lighthill  acoustic  analogy  approach. 

The  analysis  of  the  acoustic  radiation  from  rotating  blades 
begins  with  the  inhomogeneous  acoustic  wave  equation  developed  by 
Lighthill  [2],  together  with  the  boundary  condition  imposed  by  the 
blades.  The  implications  of  employing  the  Lighthill  acoustic  analogy 
approach  are  as  follows:  (1)  the  observation  point  of  the  acoustic 
radiation  is  far  from  the  perturbed  region  so  that  the  propagation  and 
source  terms  are  separable  in  the  governing  equation  for  wave  propa¬ 
gation,  and  (2)  the  medium  at  the  observation  point  is  assumed  to  be 
at  rest.  The  quiescent  medium  and  the  assumption  that  the  acoustic 
process  is  isentropic  simplify  computing  the  acoustic  intensity. 

The  Lighthill  acoustic  analogy  approach  leads  to  an  integral 
equation,  with  Green's  function,  for  3n  infinite  medium  where  the 
integration  is  done  over  three  spatial  dimensions  and  one  temporal 
dimension.  Here,  the  time  is  regarded  as  a  parameter;  i.e.,  the  so 


13 


called  recarded  Cine  concepc  is  employed.  However,  because  of  Che 
exiscence  of  che  boundary  surface  imposed  by  che  blades,  Green's 
funccion  of  che  acouscic  wave  equacion  in  an  infinice  medium  connoc  be 
used  wichouc  modifying  che  boundary  conaicion  wich  a  Heaviside 
funccion.  The  modificacion  creaces  new  source  earns  such  as  monopole, 
dipole,  and  quadrupole  sources.  However,  che  acouscic  wave  equacion 
modified  by  che  Heaviside  funccion  yields  che  same  Green's  funccion  as 
in  che  wave  equacion  of  an  infinice  medium,  so  chac  Che  wave  equacion 
wich  boundary  condicions  can  be  solved  wichouc  seeking  a  new  Green's 
funccion. 

2.2.  Incearal  Eouacion  for  Acouscic  Pressure 

Farassac  [34]  showed  chac  che  surface  on  a  moving  body  can  be 
creaced  as  a  poinc  dipole  source  whose  sCrengch  is  equal  Co  che  cocal 
force  on  che  vircual  radiacion  surface,  when  che  following  relacions 
are  sacisfied: 


r  >>  —  and 
P  a 


(2.1) 


where  t  ,  ,  and  L  are  Che  characceriscic  cime  scales  of  che 

P  Mr 

pressure,  che  Mach  number  of  Che  blade  rocacional  speed,  and  Che 
lengch  scale  of  che  vircual  surface. 

Thus,  che  acouscic  wave  equacion  Co  be  solved  for  rocacing 
blades  is 


1 

□  " 


(2.2) 


where 
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B  M 

E  E  i(x )•:(>'. 
m=l  n=l 


1R 

x -  r  )  :  (Z)  :  (R-R  )~- 
m  n  R 


(2.3) 


Fig.  1  shows  Che  configuration  of  the  relation  between  the  rotating 
dipole  and  the  observer's  point,  where  spherical  coordinates  and 
cylindrical  coordinates  are  used  for  the  observer's  point  and  the 
rocacing  dipole,  respectively. 

With  the  aid  of  the  fundamental  solucion  of  Eq .  (2.2),  i.e., 
Green's  function  of  the  acoustic  wave  equation  in  an  infinite  medium, 
one  obtains 


(p(r,3,?,c)  =7 T  ;7j7  (V?)*:(x-t-4>kdRd:dZdx 

Voi.  0  a  3 


where 


2  ->  2  1  / 

D  =  (r  +  R“  -  2rR  sir.:  cos(x-))  -  2r2  cost  +  Z  ) 


(2.4) 


(2.5) 


=  r  -  R  sin:  cos(x-3)  -  rZ  cos: 


(2.6) 


The  integration  in  Eq.  (2.4)  is  within  the  stationary  reference  frame 
fixed  to  the  center  of  the  rotor  plane.  3v  introducing  a  reference 
frame  rotating  with  che  speed  of  the  rotor,  i.e.,  3  =  .Ex+>  ,  one 
obtains  the  following  equation. 


P(r,5,:,c)  = 


I'll 


Vol . 


■f) 5(x-t+^)^dRd;dZdx 

el  u 


where 


3  =  Ex  + 


(2.7) 

(2.3) 


Using  the  relation  7*(gf)  =  g? • f  +  f*7g,  where  g  and  f  are  scalar  and 


Fig.  L  Observer/Source  Coordinates 


16 


vector  functions,  respectively,  one  obtains. 


5(x-t-*-^-)  5(t-t+^-)  5(x-t+— •) 

f)  =  - 5^-  7-f  +  £.7 - r^- 


(2.9) 


Utilizing  Eq.  (2.9)  in  Eq.  (2.7),  and  using  the  fact  that 


5  (t-CH — )  o(T-tH — ) 

I  7.(, - -a-f  )dv-J  - — 


f*ndS  =  0  ,  (2.10) 


Vol . 


Surface 


because  there  always  exists  a  surface  where  f  =  0,  one  obtains. 


P(r ,  9  ,  ? ,  t )  = 


-  ^  5(x-t+^) 

17  ////  f‘70( - ^~)RdRd>d2d: 


Vol. 


D 


(2.11) 


Taking  the  gradient  of  5 (r-t+D/a) /D ,  then  Eq .  (2.11)  become? 


00 

P(r,3,p,t)  =  -  ±Jff  / (f --qD){ 


5’(T-t+^)  5(T-t-3) 

a  a 


Vol.  -» 


aD 


,2  ' 


RdxdRdMZ  , 


(2.12) 


where  the  interchange  between  dx  and  RdRdidZ  is  possible  because  the 
integration  is  within  a  rotating  frame. 


(V  D)  =  (R  -  r  sino  cos(r-v))/D  , 
0  iv. 


(2.13) 


(V  D) =  (-r  cos  o  +  Z)/D  , 


(2.14) 


("  D)  .  =  -(r  sino  sin(i-'))/D  , 


(2.15) 


and 
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5' 


(r-t+^) 

a 


d{  o  (t-c+^)  ;• 
_ a 

d(:-cA 

a 


(2.16) 


Eq.  (2.12)  is  a  general  expression  for  acousric  pressure  due  ro 
a  distributed  force,  and  it  includes  the  oressure  i  n  the  near  and  far 
fields.  In  che  following  section,  the  tern  for  the  far  field  pressure 
is  extracted  from  Eq.  (2.12). 


2.3.  Far  Field  Acoustic  Pressure 

Utilizing  Eq.  (2.3)  in  Eq .  (2.12),  and  integrating  with  respect 
to  dR,  dZ,  and  d7, 


3  M 


f(r,:,t)  =  -  —  ;  I  :  (l  *7D  )2R  ■■ 


■  I  r  mil, 

;  (r-f - ) 


1  __ i  '  an  :  an'  "n  aD 

n=l  mn 


r  /  ,  an, 

6(t-cT' 


:-dt  , 


(2.17) 


where 


D  -  r  -  R  sino  cos(.2t+J  )  , 
mn  n  m 


(2.13) 


(7  D 
o  mn 


) 


r  sino  sin(f!T-K>  ) 
_ m 

D 

mn 


(2.19) 


and 


(7  D  )_ 
o  mn  Z 


(7  D  )„  = 
o  mn  R 


r  cos ; 
D 

mn 


R  -  r  sini  cos(r.t+i>  ) 
n  n 


(2.20) 


(2.21) 


In  the  above  equation,  9  =  0  is  assumed,  because  the  transient 
pressure  is  not  considered  here,  as  seen  from  the  integration  range 
of  t  in  Eq.  (2.17).  Mow,  integrating  Eq.  (2.17)  with 


respect  to  t 


and  imposing  the  condition  for  the  acoustic  far  field,  2rr  >> 
wavelength  (Appendix  I) ,  then  one  obtains 
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Utilizing  Eq.  (2.13)  in  Eq .  (2.22),  one  obtains 


iP(r,  i,  t) 
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3  Z  -  r 

where  1'  and  1  are  the  3  and  Z  comDonent  of  1  (t),  and  1  =  0  is 

mn  mn  •  mn  mn 

assumed  on  a  blade. 

Furthermore,  imposing  the  condition  for  the  geometric  far  field, 

r  >>  R^  for  all  n  =  1,2,...,M,  and  because  the  time  history  of  the 

pressure  is  not  considered  in  the  present  study,  one  can  neglect  the 

terms  with  R  /r  in  Eq .  (2.23),  yielding 
n 


(2.23) 
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D  /a  . 
mn 


3  Z 

Introducing  2,  (x)  «  -  2,  (t)  sinS  and  l  (t)  =  l  (x)  cos3  . 

mn  ran  ran  ran  mn  ran 

where  l  (x)  is  assumed  to  be  normal  to  the  chord  line  of  a  thin 
mn 

airfoil,  and  2  (t)  is  obtained  bv  Sears'  function  (see  Fie.  2), 

mn 


Eq.  (2.24)  becomes. 
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Eq .  (2.25)  is  an  extension  of  the  result  of  Lowson  [7]  and 
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Morfey,  et  al.  [14]  which  deals  wich  the  acoustic  radiation  from  a 

moving  point  force.  As  will  be  shown  below,  this  relates  the  acoustic 

intensity  spectral  density  function  at  the  observation  point  x  with 

the  power  spectral  density  function  of  the  turbulence  experienced  by 

rotating  blades  calculated  by  an  aero/hydrodynamic  function. 

As  shown  in  Eq .  (2.25),  the  variation  of  the  distance  between 

the  observer  and  the  blade  on  a  time  axis  causes  a  phase  modulation  in 

1  (t)  because  r  =  t  -  D  /a.  Furthermore,  the  variation  of  the  Mach 

mn  mn 

number  of  the  blade  velocity  in  the  direction  of  the  observer,  i.e., 

(R  2/a)  sinosin(._'.T+;>  ),  behaves  as  an  amplitude  modulating  signal  of 
n  m 

i  ( t )  and  11  ( T ) / 3 1 . 

mn'  mn 

The  Fourier  coefficients  of  the  modulation  terms  are  derived  in 
Appendix  II.  The  following  is  the  final  result  taken  from  Appendix 
II: 
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where  che  approximation  is  made  for  J^tz)  =  (1/p  •' )  (z/2)  p  in  Ref.  [28] 
Hence,  the  higher  order  terms  can  be  neglected  for  small  (R^O/a) sino . 


2.4.  Acoustic  Intensitv 


One  can  rewrite  Eq.  (2.25)  as  follows: 


BN00 


3Z 
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and 


r(]S|+l)r( jkj+1)  ’ 


P  =  s  +  k  ,  (2. 33) 
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Utilizing  the  following  definition  of  the  autocorrelation  func¬ 
tion  of  f(r,  0,0  , 


^  (r.o.t.a)  =  Et  ?*(r,:,t)P(r,:,t+-j);  , 
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(2.35) 
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one  obtains 
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The  instantaneous  acoustic  intensity  is  given  by  the  following 
relation: 


T(r, s,e,t) 


(2.38) 


where  the  acoustic  particle  velocity  u(r,;,9,t)  is  given  by  the 
gradient  of  P(r,$,c)  as  follows: 


u(r 


3,t) 


rjp  y  Jp_  -r  _2L  T  1 
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(2.39) 


As  r  increases,  u(r,b,3,t)  is  dominated  by  which  is  the  r- 
component  of  u.  Thus,  Eq .  (2.39)  becomes 


u(r,a,3,t)  - 


}P  - 


-  r 


(2.40) 


Hence,  the  acoustic  intensity  at  a  large  distance  r  is  given  by 
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I(r,?,6,t)  =  I  (r,J,0,c)  =  \ 
r  - 


(2.41) 


Now,  caking  the  expectation  of  the  above  intensity,  one  obtains 


Ell  (r.j.e.t)}-  k  Mr"  £i:?(r,J,5.t)P*(r,b,e,t)})  .  (2.42) 


Since  the  pressure  in  Eq .  (2.25)  is  real,  Eq.  (2.42)  becomes 


Ell  (r,*,9,c)}  =  ^  H(T( r.b.e.c) 


(2.43) 


Eq.  (2.43)  is  equivalent  to  J5  ^  (r , 5 , C , i) /ca  when  a  =  0.  Consider 
the  following  pair  of  Fourier  integrals: 


S*  (r,d,t,w)  =  j  &  (r,i,c,w)  — — - 


ca 


and 


$^(r,d,t,:a)  =  ~  /  Sj,?(r,  i,t,w)e^“JC‘da  . 


(2.44) 


(2.45) 


Thus,  from  Eq .  (2.36),  one  obtains 


1  ” 

Ellr(r,i,t)}  =  —  J  S?p(r,5,t,w)dw 


(2.46) 


accordingly,  S-^j,  ( r ,  ;  ,  t ,  w)  can  be  regarded  as  3n  instantaneous 
intensity  spectral  density  function. 

Hence,  by  taking  the  a-Fcurier  transform  of  Eq .  (2.36),  one 
obtains  the  following  instantaneous  acoustic  intensity  spectral 


densitv  function: 
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is  stationary,  Eq.  (2.47)  becomes. 
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Assuming  SVf (r , ? , t , j)  is  continuous  and  finite  with  respect  to  t,  and 
taking  a  time  average  of  Eq .  (2.49)  over  2t/.~. ,  one  obtains 
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CHAPTER  III 


3LADE  RESPONSE  TO  FLOW  FLUCTUATION 


3 . L .  Introduction 

In  Chapter  II,  the  acoustic  intensity  spectral  density  function 
was  derived  under  the  conditions  of  compactness,  i.e.,  the  virtual 
source  size  of  the  rotating  blades  is  equivalent  to  the  physical  size 
of  the  blades,  and  the  variation  of  the  radiation  time  over  the  blade 
is  very  small  compared  with  the  time  scale  of  the  surface  pressure  of 
the  blade  so  that  the  instantaneous  radiations  from  various  parts  of 
the  blade  reach  the  observer  simultaneously.  The  acoustic  intensity 
spectral  density  function  involves  the  computation  of  the  power 
spectrum  of  the  lift  on  a  blade.  The  present  chapter  considers  the 
derivation  of  the  impulse  response  function  of  a  blade  based  on  a  two- 
dimensional,  thin  airfoil  theory,  and  the  power  spectral  density 
function  of  the  lift  by  employing  the  correlation  function  of  the  flow 
fluctuations . 

With  the  aid  of  the  result  by  Jackson,  et  al.  [24],  the  three- 
dimensional  interaction  between  olaaes  and  flow  fluctuations  can  be 
approximated  by  a  two-dimensional  strip  theory  using  Sears'  function. 
The  lift  on  each  scrip  is  obtained  in  order  to  get  the  impulse  re¬ 
sponse  function  of  each  strip.  The  two-dimensional  strip  approach 
imposes  the  restriction  that  the  acoustic  radiation  due  to  parallel 
gusts  (the  spar.wise  wavenumber  is  zero)  is  dominant,  as  investigated 
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by  Araiet  [27].  Utilization  of  Sears'  function  imposes  the  restriction 
that  the  mean  angle  of  attack  is  so  small  that  the  effect  of  the  flow 
fluctuation  parallel  to  the  chord  line  of  an  airfoil  can  be  neglected 
as  a  second  order  effect.  Hence,  our  concern  is  with  the  response  of 
a  thin  airfoil  with  zero  mean  angle  of  attack  to  an  upwash  velocity 
fluctuation,  where  the  airfoil  is  assumed  to  be  lightly  loaded. 
Furthermore,  since  it  is  assumed  that  the  frozen  flow  fluctuation 
pattern  is  convected,  there  is  a  one-to-one  correspondence  between 
frequency  and  wavenumber.  In  addition,  the  blade-to-blade  interaction 
is  not  included  in  the  impulse  response  function  of  a  blade. 

3.2.  Transfer  Function  of  31ade 

In  general,  a  flow  velocity  fluctuation  with  convection  velocity 
U^  can  be  expressed  by  the  following,  four-dimensional,  Fourier  trans¬ 
form  pair: 


V (y , t )  =  — —  r  |f/jA(k  ,k,,k  w)eJ'k'y  (klLl‘Sj)T'dudk.dk,dk^  , 
(2-0*  1  “  J  1-3 

(3.1) 


and 


A(k1,k7,k3,w)  =  //J7v( y,r)e  J  *  -'’+<'^1"1+,J-  '  "dtdy1dy2dy3  . 


(3.2) 


When  the  flow  has  a  property  such  that  ~  <<  k^U  in  Eq .  (3.1), 
Eq.  (3.1)  becomes 


V(v,t)  =  — —  fr//  A(k, ,  k, ,  k. ,  j  )  e j  ( 1 T )  du.dk,  dk, 
(2-)*  *  '  J  1  “ 


dk3  . 
(3.3) 
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Since  :  A(k  ,k-,k  ,,w)du  =  A(k,  ,k  ,k-,r),  when  t  =  0,  Ea .  (3.3) 

'  “a>  I  Z  j  X  -  J 

becomes 


V ( y , r )  = — —  t  ///  A(k, ,k„,k„,r=0)e 
(2-)J  —  *  "  J 


j(k.y+k1-Ji:)dk^k?dkj  , 

(3.4) 


where  A(k.  ,k?  ,k  ,  t)  =  j  A(k^,k,  .k^  ,  j)e-' ** ‘dm.  Hence,  k.U  >>  j 

is  equivalent  to  a  frozen  velocity  fluctuation  pattern,  since 
A(k^,k2»k_,T“0)  is  the  wavenumber  spectrum  of  the  velocity  fluctuation 
at  an  initial  time  r=0.  In  other  words,  the  frozen  velocity  fluctu¬ 
ation  pattern  assumption  is  valid  when  the  typical  time  scale  of  the 
fluctuation  with  no  convection  is  much  larger  than  that  of  the  frozen 
fluctuation  with  convection. 

Due  to  the  frozen  fluctuation  pattern  assumption,  the  Fourier 

i  1c  U  t 

component  of  the  temporal  variation,  i.e.,  eJ  11’,  can  be  related  to 

j  1^  •  » * 

that  of  a  spatial  variation,  i.e.,  e  '  ,  through  k^ .  Thus,  a  one-to- 
one  correspondence  exists  between  the  frequency  and  wavenumber  domain, 
i.e.,  fc  and  k^. 

In  addition,  the  above  restriction,  j  <<  k^U^,  can  be  relaxed 
because  in  the  computation  of  the  lift  on  a  blade  it  needs  only  be 
assumed  that  the  temporal  variation  is  negligible  during  the  convec¬ 
tion  time  for  the  fluctuation  to  pass  through  the  blade. 

Assuming  linearized,  small  perturbation  analysis  of  the  thin 
airfoil  response  to  a  frozen  velocity  fluctuation  pattern,  with  a 
convection  velocity  in  the  y ^  direction,  3S  shown  in  Fig.  3,  the 

response  can  be  expressed  by  the  sum  of  the  response  to  each  Fourier 

T  f  {C  *  ^  ^  ^ 

component  A(k^,k0,k^)eJ  ■’  11',  where  A(k^,k,,k^)  implies 
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ACk-.k^.k-.T  =  0).  The  lift  on  a  blade  with  a  spanwise  length  b  and 
I  Z  j 

a  chordwise  length  c  can  be  given  by  the  following  integration, 


®  c/2  b/2+o 

VT>  -  ///  /  /  / 

-ao  -c/2  -b/2  -o 


G3. (k,y)A. (k)eJ (k'y+klLlT)5(y3)dydk  , 

(3.5) 


where  G^  (k,y)  is  an  aero/hydrodynarnic  function  in  tensor  notation. 

Neglecting  the  lift  due  to  the  velocity  fluctuations  parallel  to 
the  blade  surface,  because  of  their  second  order  effect,  and  assuming 
incompressible  flow  so  that  k  and  A  are  orthogonal  each  other,  one 
obtains 


°°  c/2  b/2+o 

VO  -//////  G 

-®  -c/2  -b/2  -o 


33  (k,v)A3  (k)  e-'  ^k  ''+klk  1 "  ^  5  (y3)dydk  . 

(3.6) 


The  integration  of  Eq.  (3.6)  over  y3  yields 
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By  assuming  G^^k^.’rf^.k^y^.y,,)  to  be  constant  over  the  blade  span, 
and  not  considering  edge  effect  (the  tip  of  the  airfoil),  Eq .  (3.7) 
becomes 
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Its  approximate  form  is 
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where  this  approximation  is  done,  not  by  taking  the  limitation  of 
Eq.  (3.12)  as  k^  increases  or  decreases,  but  rather  by  observing  that 
Eq.  (3.12)  can  be  approximated  to  within  a  few  percent  over  most  of 
its  range.  This  approximation  of  the  amplitude  was  done  by  Liepmann 
[35],  and  approximation  of  phase  was  done  by  Giesing,  et  al.  [36], 

The  integration  with  respect  to  k^  in  Eq .  (3.11)  yields  the 
contribution  due  to  the  upwash  at  y  ^  =  0,  so  that  A^(k^)  can  be 
obtained  by  measuring  the  upwash  of  the  airfoil  at  =  0,  where 

rco  4  If  jr 

A3^kl^  =  J  -a.  A3(k^,k3)eJ  3^3  dk^  at  y^  =  0.  By  using  the  notation 
A^(k^),  Eq.  (3.11)  becomes 


i3(r)  =  2bH(k1)A3(k1)ejklUlT  . 


(3.14) 


In  the  above  equation,  H(k^)  can  be  regarded  as  a  transfer  func- 
tion  for  the  response  of  a  thin  airfoil  whose  input  is  A^(k^)eJ  IT, 
Furthermore,  since  k^U^  can  be  regarded  as  an  angular  frequency,  H(k^) 
and  A- (k,  Je^l^'l  ‘  can  be  written  as  H(_/'L\)  and  A,  (w/U,  )  eJ  JT  , 
respectively,  where  w  =  k^U^.  The  implication  of  Eq .  (3.14)  is  that 
there  is  a  one-to-one  correspondence  between  a  wavenumber  and  an 


angular  frequency,  that  Che  airfoil  can  be  created  in  cine  and  fre¬ 
quency  domains  instead  of  space  and  wavenumber  domains,  and  that 
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Sears'  function  can  be  regarded  3s  the  t-Fourier  transform  of  the 
impulse  response  function  of  the  airfoil. 

Accordingly,  Eq.  (3.14)  yields  the  following  convolution 
integral, 

CD 

*3(t)  =  /  h(T-a)u(a)da*2b  ,  (3.15) 

—  CO 

where  u(a)  is  the  upwash  measured  on  a  time  axis. 

Finally,  since  i^( r)  is  obcained  by  integration  over  the  entire 
airfoil  or  the  strip  of  the  airfoil  so  chat  the  airfoil  or  strip  can 
be  replaced  with  a  dipole,  the  upwash  on  the  integration  area  has  Co 
be  well  correlated  so  that  the  dipole  can  represent  the  lift  on  the 
integration  area.  As  an  important  result,  since  the  lift  is  obtained 
in  a  time  domain  and  in  a  frequency  domain,  Eq.  (3.14)  becomes  insen¬ 
sitive  to  the  geometrical  orientation  and/or  position  of  the  airfoil 
explicitly,  although  the  upwash  includes  the  geometrical  information 
implicitly . 

3.3.  Correlation  Function  of  Lift  and  its  Time  Derivative 

Adopting  the  model  discussed  in  Section  3.2.,  one  can  obtain  the 
th  _  ch 

following  lift  at  the  n  strip  or  the  m  blade. 


l  (t)  =  J  h  (r— ,)u  (-)dv.‘.X  , 

jnn  *  '"n  -in  n 


nn  rzn 


(3.16) 


where  Eq .  (3.15)  is  utilized,  and  u  ft)  is  the  upwash  ac  the 


ch 
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strip  of  the  mC^  blade.  Hence,  the  correlation  function  between 
£^("0  and  2.^(1),  where  k  and  1  denote  the  kCh  blade  and  1th  strip, 
respectively,  is  given  by 


^mnki(T’a)  =  EUmn(T)iki(T+a)} 


(3.17) 


Utilizing  Eq.  (3.11)  in  Eq.  (3.17),  one  obtains 


U 

mnki 


(T>a)  =  ^  hmn(nl)hkZ(n2)CkZ(T’a-n2+nl)drildn2iRnAR 


l  ' 
(3.13) 


By  taking  the  ct-Fourier  transform  of  Eq.  (3.18),  one  obtains 

Cci<T’“)AR»ARz  ■  <3-19) 

The  correlation  function  and  the  spectrum  of  l  (x)  and 

mn 

(T)  as  well  as  those  of  3Z  (t)/3x  and  31,  „(t)/9t  are  needed  to 
k*.  mn  kx. 

obtain  the  autocorrelation  function  of  the  acoustic  pressure  at  the 
observer  point.  By  taking  the  time  derivative  of  Eq.  (3.16),  one 
obtains 


31  (t)  «  , 

m  /  V\n(^''»u  (*)dn- 

jT  mn  an 

—  CO 


(3.20) 


Eq .  (3.20)  shows  that  the  time  derivative  ot  1  (t)  can  be  regarded  as 

mn 

the  output  of  the  system  whose  impulse  response  function  is  Vn  (t)/'t. 

mn 

Therefore,  one  obtains  the  relieving  by  the  sane  method  as  used  in 
obtaining  Eq.  (3.19): 
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and 


where  i 


zz'  * 

smnki(-.w)  =  H  (m)ju,H,  .(j)SUU  . (t,u)AR  *R.  , 
ran  J  ttx.  ranK.A,  nl 


(3.21) 


0.22) 


SmnkZ(T-j)  =  “2h*  (ui)H,  ,U)SUU  .(T,w)iR  iR,  , 
mnK-x*  ran  <*,  ranK4  ~  ?  * 


n  1 


(3.23) 


'  denotes  the  time  derivative  of  the  lift 


CHAPTER  IV 


POWER  SPECTRUM  OF  TURBULENCE  INTO  ROTOR 

4.1.  Introduction 

As  shown  in  Chapter  II  and  III,  the  power  spectrum  of  the  inflow 
to  a  rotor  is  necessary  in  order  to  obtain  the  acoustic  intensity 
speccral  density  function  of  the  noise  radiation  from  the  rotor.  The 
present  chapter  deals  with  the  power  spectrum  density  function  for  the 
upwash  at  rotating  blades  with  and  without  inlet  struts. 

For  turbulent  flows  with  no  struts,  two  cases  are  to  be  consid¬ 
ered:  (1)  the  turbulence  is  homogeneous,  isotropic,  and  frozen,  and 

(2)  Che  turbulence  in  (1)  is  elongated  on  its  course  to  the  rotor 
plane.  The  second  case  is  pursued  to  obtain  results  for  practical 
situations  where  the  eddy  to  a  rotor  is  elongated,  i.e.,  anisotropic. 

For  turbulent  flows  with  struts,  it  is  assumed  that:  (1)  the 
turbulence  in  the  non-disturbed  inflow  and  wakes  behind  the  inlet 
struts  can  be  treated  as  in  the  case  with  no  struts,  (2)  the  wake 
turbulence  is  statistically  independent  of  the  non-disturbed  inflow 
turbulence,  (3)  the  turbulence  in  each  wake  is  statistically  independ¬ 
ent  of  that  in  ocher  wakes,  (4)  the  random  variables,  such  as  the 
center  position  and  width  of  each  wake,  are  governed  bv  uniform  dis¬ 
tribution  functions  where  these  random  variables  are  statistically 
independent  of  each  other,  and  i'S'i  no  overlap  occurs  among  wakes. 

Furthermore,  in  the  situation  with  inlet  struts,  the  mean  flow 
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velocity  defect  behind  vanes  and  the  stator-rotor  interaction  due  to 
potential  fields,  such  as  were  investigated  by  Kemp,  et  al.  [37],  can 
be  regarded  as  a  sceadv  non-uniform  flow.  A  deterministic  function 
is  introduced  in  order  to  take  this  into  account. 

4.2.  Isotropic  Turbulence  with  Mo  Struts 

In  statistically  homogeneous  and  stationary  turbulence,  the 
correlation  function  of  the  turbulence  velocities,  u^  (Pj.  •  and 
Uj^2,C2^  where  =  *■»  3,  is  given  by 

*  -►  -*■  >wu  • u  ‘  - 

E{ui(o1,t1)Uj  (p,,t2)  ;  =(n.1^  J(p,a)  ,  (4.1) 

where  s  =  ancl  s  =  “  C2'  Furt;hermore ,  by  assuming  frozen, 

isotropic  turbulence,  Eq .  (4.1)  becomes 

E{u*(p1,t1)u.(p2,t:);  =  ^jUj(|p;)  .  (4.2) 


Under  the  above  assumptions,  i.e.,  that  turbulence  is  homogeneous, 
frozen,  and  isotropic,  von  Karman,  et  al.  [38]  derived  the  following 
correlation  function: 

.  u 


i.l  =  1.3,3  ,  (4.3) 


2  r  f  (  o  -  i  -  t  C  o  ) 


+  s(! 


where  :(  ;  )  is  a  longitudinal  correlation  coefficient  function 


defined  bv 
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u“ 


(4.4) 


g(jp|)  is  a  transverse  correlation  coefficient  function  defined  by 


\  1 

_  1  2 

g(  |o  j )  =  - — - 


(4.5) 


2 

and  u~  is  the  mean  square  value  of  three  components  of  the  velocity, 
Ul’  U2’  an<*  u3’  such  c^ac 


2 

u 


(4.6) 


throughout  the  homogeneous  and  isotropic  velocity  field. 

Furthermore,  in  an  incompressible  fluid  medium  the  two  corre¬ 
lation  coefficients,  g(|pj)  and  f(jpj),  have  the  following  relation 
as  shown  in  Ref.  [39]: 


g(|o|)  =  f(|o|)  + 


if(  0‘) 

i|o, 


(4.7) 


Now  consider  the  correlation  function  of  the  upwash  at  rotating 
blade  under  the  following  assumptions:  (1)  an  incompressible  fluid 
medium,  (2)  frozen  turbulence  pattern  convected  with  an  axial  flow 
velocity  l'^,  and  (3)  homogeneous  and  isotropic  turbulence,  see  Fig.  4. 
Under  these  assumptions,  Lane  [33]  derived  the  correlation  function 
of  the  upwash  at  rotating  blades  provided  that  the  longitudinal 
correlation  coefficient  f(|t.)  is  assumed  to  be 
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f( 


exp  i  - 


(4.3) 


and  utilizing  Eq.  (4.3)  and  Eq .  (4.7). 

The  upwash  correlation  function  derived  by  Lane  [33]  is  as 
follows : 


where 


and 


iu  ,  x  2  \T.,  ,  ,,  o(a),  ,  F (et)  -  ,  !o 

mk,(a)  =  u  {E(a)  (1  -  — )  +  1-^yr-;exp( 

l  l 


A 


-)  , 


s  2,  ,  ,  „  2tt (k-ra)  N  .  2, 

E(ct)  =  cos  3  +  cos(.ia  +  - - - )sm  3  , 

D 


2  2  °  ?  2Tfk-m) 

F(a)  =  (U,a)  cos  3  +  R  R._,sin~3  sin  (La  +  - — — )  , 

■i  nin  fCx>  d 


.2  .  „2  ,  _ 2 


(a)  =  (U,a)“  4-  R*-  +  RT  -  2R  R,  ,  cos (La 

3  mn  Tel  an  .<1 


.a  + 


2^ (k-a) 
B 


(4.9) 

(4.10) 

(4.11) 

.(4.12) 


The  stagger  angle  3  is  aoproxiaatea  bv  3  constant,  3,  for  all  a 
mn 

and  n.  In  other  words,  the  variation  of  the  blade  twist  is  assuaed 
to  be  small. 

Eq.  (4.9)  can  be  written  in  terms  of  non-diaensional 
variables,  so  that 


o.(a)  F,(a) 


dk*00  =  u2  ^(3)(1  "  JT~'  +  27T(T7'-ex?(-  “i(^>>  (i-13) 


where 


2/  ,  rU3^~  f  an 

3i(3)  -  LT~  +  —  - 

v  i  .  ^  i  - 


R  2  ,-R.  .'  2  R  R.  . 

.<  <-  ,  an 


j.  ' 


'i  -  -  . 2’ (k-a) s  , ,  ,  ,  s 
cost. -oh - — — )  ,  (4.1-.) 


and 


I U3'-)  “  2  RanAki  2 

F,  (a)  =  ,  —  ■  cos  i  +  - i - sir. 

1  i -i >  i: 


.  .2,.  ,  2 a (k-a) ,  ,,  , r, 

-  sm  i„a  +  - r - )  .  (4.15) 
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Now  consider  Che  case  chat  R^  =  R,  .  =  R  for  all  m,  n,  k,  and 
l.  Eq.  (4.13),  Eq.  (4.14),  and  Eq .  (4.15)  change  as  follows: 


ruu,  .  2,  °l(:t)  .  Fl(')  ,  -o,  (a) 

$mk(a)  =  u  ,-EU;(1  ’  2  }  +  1 


where 


and 


0i(o>  iv  ■  'I'., 


.  2” (k-m) , , 
cos  (.:a  +  - - - )  ;  , 


rU3a'2  2.  ,  TrI2  .  2.  .  2,„  2*(k-mK 

sin  :  sin  ('-.a  +  - — - ) 


(a)  -  cos  3  +  ^ 


B 


(4.16) 


(4.17) 


(4.18) 


4.3.  Non-isotrooic  Turbulence  wich  No  Struts 

It  was  shown  by  Ribner,  et  al.  [40],  based  on  the  continuity 
and  conservation  of  angular  momentum  of  fluid  elements,  how  a 
contracting  flow  changes  the  incoming  turbulence  spectrum  and  its 
intensity  and  length  scales.  This  type  of  contracting  stream  occurs 
in  the  flow  into  a  stationary  rotor.  However,  in  attempting  to  apply 
this  result,  the  spectrum  obtained  in  a  wavenumber  domain  is  very 
difficult  to  convert  to  one  in  a  frequency  domain.  This  difficulty 
arises  because  of  the  rotation  of  the  blades  and  because  the  frequency 
is  a  function  of  the  angle  between  the  wavenumber  vector  and  the  blade 
stagger  angle  orientation. 

Furthermore,  Chandrasekhar  [41]  shewed  that  the  general  form  of 
the  velocity  correlation  tensor  in  an  axis;— etr ic  turbulence  can  be 
expressed  by  two  arbitrary  scalar  functions.  However,  the  task  of 
obtaining  the  arbitrary  scalar  functions  constituting  the  correlation 
tensor  is  difficult  due  to  the  fact  that  the  power  spectral  densitv 
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function  of  the  cot. "elation  function  must  satisfy  the  non-negative 
condition.  For  example,  as  pointed  out  by  Kerschen,  et  al.  [42],  the 
restrictive  condition  is  noc  satisfied  even  by  selecting  a  one¬ 
dimensional  correlation  function,  such  as  an  exponential  function. 

Since  the  treatment  of  the  isotropic  turbulence  has  been  done 
kinematically  and  geometrically  in  Section  4.2.,  that  treatment  is 
extended,  according  to  Ribner,  et  al.  [40],  to  the  turbulence  ingested 
by  a  rotor.  The  inflow  turbulence  scale  is  assumed  to  be  modified  in 
the  following  manner: 


^ij('YXl,YX2’3x3) 


(4.19) 


where  the  superscript  "in",  5,  and  v  denotes  the  inflow  isotropic 
Curbulence,  the  turbulence  elongation  rate,  and  the  turbulence  con¬ 
traction  race,  respectively.  The  above  assumption  asserts  that  the 
statistics  governing  the  inflow  turbulence  are  not  changed  through  the 
flow  contraction  but  that  the  shape  of  the  domain  of  Che  statistics  is 
distorted.  In  other  words,  the  isocorrelation  contours  are  distorted 
through  the  flow  contraction. 

The  correlation  function  of  the  upwash  between  two  points  (m,n) 
and  (k ,1)  in  a  contracted  flow  field  is  given  by 


:,(i) 


(a)  F,,  (a)  •)  (  \ 

- 1  4.  — o -  e 

1  2c0(a);e 


where 


_  ,  .  f'V’12  2.  Rnn ^k.-.  .2.  .2  2-fk-m). 

r,(i)  =  ITT-  COS  O  +  — , — —  ain  c  sm  (..cH - - - > 


^•'•i  1 


V  ~  \ 


B 


?  2  "■  f  ^ 

E?(a)  =  cos  ~i  +  cos  ( .la  +  — ; — )sin~£  , 


(4.20) 


(4.21) 

(4.22) 
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and 


P^(ei)  = 


U,an*2  rR  'i  2  ,R, 
3  !  (  mn i  >t 


+  !-55L  + 

15  A,  J  ivA,-  j 


jar 


^mn^ki  .  2-rCk-m), 

cos  (..»  +  - - -  )  . 


■y  -> 

*"AZ 


B 


(4.23) 


When  R  =  R,  =  R  for  all  m,  n,  k,  and  l,  Che  above  equations  become 
mn  k 


Ckz^  "  <«ln>-i=2(»><i  - 


,(a)  F.u) 


where 


F2(3)  = 


(  2l 


2o2(a) 


-o  (a) 

e  2  ,  (4.24) 


V'l2  2 ,  .  (  R  l2  .  2,  .  2„.  ,  2t  (k-m)  s  ,,  ,c, 
cos  3  +  j— — :  sin  ssm  (--a  +  r - )  ,  (A. 25) 


\ 


iJ 


IvA^j 
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x  2,  ,  ,  2*rfk-m)N  .  2, 

E~(a)  =  cos  3  +  cos  (.-a  +  - - - )sm  3  , 

4  D 


(4.26) 


and 


2  fU3al 2  (  r  i ' .  ?-rfk-m)  , 

: 2 (a)  =  [JTi  +  2(7T7,  a  *  C0S(:-3  +  — 3  );  •  (4>27) 


4.4.  Power  Snectrum  of  Uowash  with  Struts 


The  upvash  to  blades  operating  downstream  of  a  set  of  struts  is 
assumed  to  be  expressed  by  the  following  equation: 


u 

mn 
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S-i 


Rect  i- 


q=0 


2-rq  ^  2m 
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?qm 


com  i  l  .  . 
"-J  ru  (t) 
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QO  S-1 

r  c 

p=-°°  q*0 
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t 

q 

— i  u 
'  mn 


(t) 


+  M 

mn 


<  t)  * 


(4.23) 


where  u1  and  are  the  turbulent  velocities  of  the  inflow  and  the 
mn  mn 
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wake,  respectively,  and  it  is  assumed  that:  (1)  no  wake  is  overlapped 
with  anocher,  (2)  the  function  Rect(i)  is  given  by 

Rect(r)  =  |l  It  j  <  1 

(0  otherwise  ,  (4.29) 


and  t  and  w  are  random  variables  for  the  position  and  width  of 
qpm  qpm  r 

the  Rect(r)  function  whose  probability  density  functions  are  given  bv 
the  following: 


and 


Pr(Cqpm)  =  fll 


t  '  <  A 
qpm'  = 

otherwise  , 


P  (w  ) 
*r  qpm 


( — - — 

j ^2""1 
(0 


<1,  <  w  1  A_ 
1  —  qpm  —  2 

otherwise  , 


(4.30) 


(4.31) 


where  A,  A^,  and  A0  are  determined  so  that  no  overlap  among  wakes 

occurs,  (3)  the  function  u^(t)  is  the  turbulent  velocity  in  the  wake 

of  the  qth  inlet  strut  interacting  with  the  nCh  segment  of  the  mth 

rotating  blade,  (4)  e{u*  (r)}  =  Z{uq  (t)j  =  0  for  all  m,  n,  k,  l ,  and 

9 i  (5)  uqn(T)'s  ^or  9  =  0. - ,S-1  are  statistically  independent  of 

each  ocher,  (6)  the  correlation  function  of  u1  (t)  and  u,1  (t)  is  given 

mn  kx.  ° 

by  Eq.  (4.20),  (7)  the  correlation  function  of  uq  (t)  and  uq(r)  is 

mn 

given  by 
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(4.32) 


where 
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F2(ct) 


ru3a  12 
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(4.33) 


E?(a)  =  cos“3  +  cos(P.a  +  2,1  ^  m-)sin^3 
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(4.35) 


(8)  the  function  M  (t)  is  a  deterministic  function  due  to  the  steady 
circulation  and  mean  velocity  defect  of  the  upstream  strut,  and 

(9)  the  variables  t^^  an^  wqpm  are  statistically  independent  for  all 
q,  p,  and  m. 

Now  consider  the  correlation  function  of  the  upwash  which  is 
given  as 
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By  caking  che  expeccation  of  Eq .  (4.36),  one  obtains 
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-xpanamg  che  tunccions  g  (t)  and  M  (~)  in  a  Fourier  series, 

~n  ’ 

Eq.  (4.39)  becones 


(4.41),  one  obtains 
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Where  Smnki(to)  is  che  Courier  transfora  of  ^^(a)  for  all  q, 
assuming  S^^Ca)  is  the  same  for  all  q  =  0,1 . S-l. 


CHAPTER  V 


ACOUSTIC  INTENSITY  SPECTRUM 
DUE  TO  INTERACTION  OF  A  ROTOR  WITH  TURBULENCE 

5.1.  Introduction 

In  Chapter  II,  the  acoustic  radiation  intensity  from  rotating 
blades  in  a  free  space  was  obtained  and  was  related  to  the  power 
spectrum  of  the  lift  of  rotating  blades  under  the  condition  chat  the 
distributed  pressure  on  each  blade  is  compact.  In  Chapter  III,  the 
impulse  response  function  of  a  blade  was  derived  so  that  the  lift  can 
be  obtained  in  a  time  domain  and  in  a  frequency  domain.  In  Chapter 
IV,  the  correlation  function  of  the  upwas'n  to  the  blades  and  its 
spectrum  were  derived. 

In  the  present  chapter,  the  results  in  Chapter  II,  III,  and  IV 
are  combined  so  that  the  acoustical  intensity  spectrum  can  be 
obtained.  In  deriving  the  explicit  expression  for  the  acoustic 
intensity,  it  is  shown  that  the  retarded  time  difference  of  the  lift 
among  blades  can  be  neglected.  The  condition  necessary  for  this 
demonstration  to  be  valid  is  more  relaxed  than  that  for  the  dipole 
domination  over  other  radiation  sources. 

5.2.  Retarded  Time  Difference 

Since  the  results  in  Chapter  IV  were  obtained  at  the  time 
associated  with  the  emission  of  the  acoustic  radiation  on  blades, 
those  results  cannot  be  substituted  into  the  results  in  Chapter  II 
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where  the  correlation  is  taken  at  the  observation  time  when  the  radi¬ 
ation  reaches  the  observer's  point.  In  other  words,  the  correlation 
function  for  the  blade  lift  S^\’^(',a)  i-s  obtained  at  two  different 

times  t  -  D  /a  and  t  +  a  -  D,  , /a. 

mn  iti 

Now  consider  under  what  condition  the  following  approximation  is 

attainable  so  that  the  retarded  time  difference  can  be  neglected 

in  the  correlation  function,  $  (t  ,a) ,  in  Eq.  (2.36), 

mnk.  a# 
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mn 


(t 
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l  (t  -  -) 
mn  a 


(5.1) 


By  Taylor  series  expansion,  one  can  get  the  following  expansion 

of  H  (t  -  D  /a), 
mn  mn  ’ 

D  di  (t  -  -) 

„  .  mn.  „  .  r.  ,  mn  a 

l  (t - )  =  l  (t - )  +  - - -  < 

mn  a  mn  a  dt 


,  r  ~  •  nsris  .  ' 

R  sini  cos-,  ..(t - ;  +  >  ; 

c  an  am' 

- ;  +  error  terms  . 


(5.2) 


Hence,  if  it  can  be  assumed  that  the  typical  time  scale  of  1  (t)  is 

mn 

very  large  compared  with  R  /a  or>  equivalently,  that  the  upper  cutoff 
frequency  of  ^(r)  is  much  lower  than  a/R  ,  then  the  approximation 
in  Eq.  (5.1)  is  attained.  This  condition  is  relatively  relaxed  when 
compared  with  the  one  for  the  dipole  domination  over  other  radiation 
sources  such  as  was  given  by  Goldstein  [6]: 

L 

a(l  -  Mr)  ’ 


T  >> 

m 


(5.3) 
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where  t  ,  L,  and  M  are  a  characteristic  time  for  the  source  fluctu- 
m  r 

ations  measured  in  the  moving  frame  with  a  convected  flow,  a  charac¬ 
teristic  source  dimension,  and  a  Mach  number  of  the  body  speed 
component  in  the  direction  of  the  observer,  respectively.  Further, 
fans  with  low  rotational  speed  and  small  radial  sice  are  more  likely 
to  satisfy  the  condition  for  the  approximation  in  Eq.  (5.1). 

Since  A  (t  -  D  /a)  is  given  bv  the  convolution  integral  in 
mn  mn 

Eq.  (3.15),  then  the  following  equation  is  obtained  so  that  the 

condition  required  for  the  approximation  in  Eq.  (5.1)  is  imposed  on 

either  h  (r)  or  u  (t) . 
mn  mn 


h  (t-n)u  (n)dn 
mn  mn 


h  ( n ) u  (-r-n)dn 
mn  mn 


(5.4) 


Thus,  one  can  obtain  the  condition  for  the  approximation  to  be  imposed 


on  h  (t)  or  u  (t);  i.e.,  the  typical  time  scale  of  h  (r)  or  u  (t) 
am  mn  mn  mn 

should  be  much  larger  than  R^sinb/a  or  equivalently,  the  upper  cutoff 

frequency  of  h  (t)  and  u  (t)  are  derived  in  the  same  fashion  as  was 
ran  mn 

done  in  Eq.  (5.2).  This  condition  is  easily  obtained  when  the 


observer  is  near  the  axis  of  the  rotor,  and  is  relaxed  when  compared 


with  the  condition  of  dipole  domination,  (5.3). 


First,  consider  the  condition  on  u  (t),  using  its  typical 


length  scale  L  and  its  convection  velocity  R  P,  where  l’,  <<  R  P  is 

'  mn  3  mn 

assumed.  According  to  the  condition,  one  obtains  the  following 


expression: 


R 

mn  a 
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L  R  siro 
mn 


(5.5) 
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Thus,  low  rotational  speed  and/or  small  ratio  of  rotor  radius  and 
turbulence  length  scale  are  required  so  as  to  obtain  the  approxima¬ 
tion.  As  mentioned  before,  the  condition,  (5.5)  is  more  likely  to  be 
satisfied  in  water  than  in  air.  When  the  observer  is  at  J  °  0,  the 
condition  is  automatically  satisfied,  i.e.,  there  is  no  retarded  time 
difference  for  the  lift  distribution  on  the  blades. 

Next,  consider  the  condition  on  h  (t).  Let  f  be  the  upper 

mn  cut 

cutoff  frequency  of  h  (t) ,  where  f  is  defined  bv  -3  dB  =  10  log 
H  mn  cut  & 

(H  ( 2tt f  )/H  (0)}.  Then,  the  condition  to  be  imposed  on  h  (t)  is 

mn  cut  mn  mn 

obtained  as  follows. 
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(5.6) 


(5.6)  shows  that  low  rotational  speed  and/or  low  aspect  ratio 

rotors  are  most  likely  to  satisfy  the  condition. 

Furthermore,  for  the  time  derivative  of  i  (t),  the  results  for 

mn 

u  (t)  and  h  (t)  are  applicable  based  on  the  following  equations: 
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(5.3) 


I 

In  the  present  study,  either  one  of  the  results  (5.5)  and  (5.6), 
or  both  of  them  mav  be  assumed  to  be  satisfied  by  the  rotor  such  that 
the  rotor  is  operated  at  low  rotational  speed,  the  flow  fluctuation 


length  scale  is  large  compared  with  the  rotor  radius,  and  the  cutoff 
frequency  satisfies  (5.6).  In  the  subsequent  derivations,  the 
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retarded  time  difference  is  neglected  in  the  spectrum  analysis  for  the 
rotor  in  a  free  space.  It  should  be  noted  that  this  assumption  is 
imposed  only  on  the  retarded  time  difference  associated  with  the  lift 
on  a  blade,  while  the  retarded  time  difference  associated  with  the 
distance  becween  the  rotor  and  the  observer  is  included  in  the  present 
study  as  was  done  in  Chapter  II. 

5.3.  Acoustic  Intensity  Spectrum  from  a  Rotor 

Operating  in  Free  Space 

By  neglecting  the  retarded  time  difference,  under  the  assumption 
discussed  in  Section  5.2.,  Eq .  (2.36)  becomes 
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(5.9) 


Then,  caking  the  a-Fourier  transform  of  Eq .  (5.9),  and  using  Eq . 
(3.19)  through  Eq.  (3.23),  one  obtains 
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Utilizing  Eq.  (5.10)  in  Eq.  (2.47),  and  replacing  t  -  r/a  with  t,  one 
obtains 
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(5.11) 


Taking  the  time  average  of  Eq.  (5.11)  over  the  period  2tt/Q,  under  the 

assumption  that  Sjr  (r,?,t,w)  is  continuous  and  finite  with  respect  to 

t,  introducing  the  typical  section  idea  such  that  H  (u)  =  H(oj)  for 

mn 

all  m  and  n,  and  replacing  a  blade  with  a  dipole,  one  obtains 
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Before  utilizing  -q .  (-*.*■3)  in  cq .  (5.12),  let  us  consider  the 
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spectrum  S^Cjj)  can  be  expressed  as 


2  U,  R  R, 

Suu(m)  =  H_  s/_l  _3L  ^ 
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(5.13) 


Mow,  utilizing  Eq.  (4.13),  Eq.  (5.13),  and  Ec .  (3.13)  in  Eq.  (5.12), 
and  non-dimensionalizing  the  result,  yields  the  following: 
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and  R^.'  -  { (R^.3)  "+U“ }  is  assumed. 


When  taking  the  turbulence  contraction  and  elongation  into 
account,  one  obtains 
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5.4.  Acouscic  Intensity  Spectrum  from  Open  Duct 
via  Free  Rotor  Model 

This  model  predicts  the  free-space  radiation  and  hence  may  be 
regarded  only  as  a  source  model.  Because  of  this  simplicity,  the 
effect  of  the  duct  is  not  involved  in  the  computation.  As  pointed  out 
by  Lansing  [43],  this  model  appears  to  overestimate  greatly  the 
radiated  noise,  and  to  have  a  completely  different  directivity  pattern 
than  his  exact  theory  claims.  However,  this  model  can  be  used  for 
diagnostic  purposes  on  fans. 

Again,  neglecting  the  retarded  time  difference,  as  was  done  in 
Section  5.2.,  and  utilizing  Eq .  (4.45)  in  Eq.  (5.11),  Eq.  (5.11) 
becomes  as  follows: 
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Now,  taking  Che  tine  average  of  £q.  (5.17)  over  2rr/72,  one  obtains 
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where  ic  is  assumed  chat  che  terms  e^,  f^  for 

n  2  n  2 

(q  -  p)  >  S  are  negligible,  and  that  Syj>  (r ,  p ,  t  ,u)  is  finite  and 

continuous  with  respect  to  t.  In  fact,  as  (R  P.sinp)  /a  decreases  the 

n 

above  assumption  becomes  more  accurate. (See  Eq.  (2.30)  and  Eq .  (2.31).) 

Now  representing  a  blade  by  a  dipole,  and  introducing  the  typi¬ 
cal  section  idea  for  the  impulse  reseponse  function  of  the  blade,  one 
obtains 
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Utilizing  Eq .  (5.15)  in  Eq.  (5.19),  and  non-dimensionalizing  the 
result,  one  obtains 
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When  there  is  no  strut,  setting  a^s)=  b(s)=  0  for  all  s,  in  Eq. 


(5.20)  and  Eq.  (5.21),  yields  an  equation  equivalent  to  Eq .  (5.14)  and 
Eq.  (5.16),  respectively. 

When  the  observer  is  at  J  =  0,  then  Eq .  (5.21)  becomes 
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CHAPTER  VI 


NUMERICAL  RESULTS 


6.1.  Introduction 

In  Chapter  V,  the  time  averaged  intensity  function  of  the 
acoustic  radiation  from  rotating  blades  was  derived  by  replacing 
blades  with  rotating  dipoles.  This  derivation  was  done  under  the 
following  turbulent  inflow  conditions:  turbulent  ingestion  with  no 
inlec  strut  wakes,  inflow  turbulence  elongation  and  contraction  with 
no  inlet  strut  wakes,  and  turbulent  ingestion  with  inlet  strut  wakes. 

For  the  present  derivation,  a  blade  was  assumed  to  be  acousti¬ 
cally  compact.  A  rotating  dipole  model  was  pursued  so  that  the  effect 
of  the  dipole  rotation  can  be  accounted  for  explicitly.  Sears'  func¬ 
tion  was  regarded  as  an  impulse  response  function  of  a  blade,  so  that 
Che  analysis  was  done  in  time  and  frequency  domains.  There  was  no 
blade-to-blade  interaction  to  be  included,  such  as  that  developed  by 
Kemp,  et  al.  [1].  Further,  the  dipole  source  was  assumed  to  be 
predominant  over  other  radiation  sources,  so  that  the  retarded  time 
difference  among  the  lift  of  blades  can  be  neglected.  It  was  assumed 
that  the  inlet  struts  do  not  affect  the  sound  propagation  and  that  the 
distance  from  che  rotor  to  the  duct  exit  is  much  shorter  than  the 
sound  wavelength  so  that  duct  effects  on  propagation  are  negligible. 
Finally,  the  turbulence  was  assumed  to  be  frozen,  homogeneous,  and 


isotropic. 
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In  this  chapter,  computations  are  done  for  the  results  in 
Chapter  V  such  as  Eq.  (5.16)  and  Eq.  (5.22)  by  changing  the  number  of 
struts  and  rotor  blades,  the  spacing  between  struts  and  a  rotor,  the 
observation  angle,  and  non-dimensional  parameters  such  as  U^/A^  and 
R  /Ai>  In  particular,  the  variation  of  these  non-dimensional  param¬ 
eters  is  made  around  an  axial  flow  velocity  of  =  30  m/sec,  a  rotor 
radius  of  R^  =  0.30  m,  a  rotor  rotational  speed  of  H  =  220  rad/sec,  an 
inflow  turbulence  length  scale  A^  =  0.15  m  and  0.017  m,  and  a  wake 
turbulence  length  scale  A  *  0.0068  m. 

6.2.  Numerical  Results 

Figs.  5,  6,  and  7  show  the  non-dimensionalized  (S  j,y(r  ,  J,  t  ,u)^> 
dependent  on  the  non-dimensional  parameters,  U^/ A^fi  and  R^/A^.  These 
figures  show  that  as  U^/A^Q  decreases,  the  level  of  the  blade-passing 
frequency  and  its  multiples  increases,  and  as  R^/A^  increases  those 
levels  also  increase.  Further,  as  R^,/ A^  increases,  the  overall  level 
increases . 

Hence,  it  is  clear  that  the  kinematic  relation  R^/A^  cannot 

explain  fully  the  spectrum  profile  without  the  dynamic  relation 

U./A.O.  The  effect  of  U-/A.0  on  the  spectrum  is  the  same  as  that 
3  i  3  i 

obtained  by  Homicz  [8].  However,  the  effect  of  R^/A^  on  the  spectrum 
for  constant  U^/A^.T  is  the  opposite  of  that  claimed  by  Mani  [15]. 

This  is  probably  due  to  the  fact  that  Mani  did  not  consider  the  vari¬ 
ation  of  U./A..7  when  he  varied  R_/A.  as  a  oarameter.  Further  investi- 
3  l  i.  x 

gation  of  this  macter  is  necessary. 

Further,  as  the  rotational  speed  increases,  the  bandwidth  and 
Che  level  at  the  blade-passing  frequencv  and  its  multiples  are 
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increased,  and  as  the  axial  flow  velocity  increases,  the  blade-passing 
frequency  and  its  multiples  appear  clearly.  The  above  phenomena  can 
be  observed  through  the  variation  of  I'y' in  Figs.  5,  6,  and  7. 

Fig.  8  shows  that  as  B  increases,  the  bandwidth  of  the  peak  at 
Che  blade-passing  frequency  and  its  multiples  increases.  In  this 
example,  as  B  is  increased  from  3  to  9,  the  width  at  the  blade-passing 
frequency  has  approximately  triples.  In  this  figure,  the  horizontal 
axis  is  w/BO.  This  tendency  is  mentioned  also  by  Homicz  [8]. 

Figs.  9,  10,  and  11  show  the  observer  angle  dependence  of  the 
non-dimensionalized  yS yj,  (r ,  p,  t ,  w)-)  .  Fig.  10  shows  that  the  domination 
of  the  blade-passing  frequency  is  reduced  as  the  observer  angle 
increases.  The  extra  wiggles  in  Fig.  10  are  due  to  the  dipole  rota¬ 
tion.  This  shows  that  the  disk  model  of  a  rotor  is  not  sufficient 
even  when  the  rotor  radius  is  about  0.3  and  its  rotational  speed  Mach 
number  is  about  0.2. 

In  the  computations  in  Figs.  9  through  11,  the  summation  in 
Eq.  (5.16)  is  done  for  the  first  five  terms  after  observing  that  the 
computed  result  from  the  first  seven  terms  is  not  different  from  that 
using  the  first  five  terms. 

Fig.  12  shows  the  influence  of  the  kinematic  contraction  and 
elongation  of  the  inflow  turbulence  on  the  noise  spectrum.  The  elon- 
g3Cion  and  contraction  produce  a  large  effect  on  the  spectrum  as 
reported  by  Hanson  [16];  although  Hanson's  elongation  and  contraction 
are  done  by  the  statistics  of  the  eddies  chopped  by  rotating  blades. 

As  shown  in  Fig.  L2,  the  elongation  and  contraction  rates,  2.0  and 
0.5,  respectively,  create  peaks  at  the  blade-passing  frequency  and  its 
multiples  even  if  there  is  no  high  level  of  peak  with  no  contraction 
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and  elongation. 

Figs.  13,  14,  15,  and  16  show  the  non-dimensionalized 

(r ,4>,t,ui)')  for  the  rotor  operating  in  strut  wakes.  In  this 

example,  there  is  no  clear  influence  of  the  wake  turbulence  on  the 

spectrum  profile.  However,  the  effect  of  the  mean  flow  deficit  due  to 

the  inlet  struts  is  pronounced  unless  the  number  of  struts  and  rotor 

blades  are  chosen  such  that  the  least  common  multiple  of  B  and  S  times 

ft  is  in  the  cutoff  range  of  Sears'  function. 

In  Figs.  14  through  16,  the  wake  position  function  Rect(t)  is 

assumed  to  be  deterministic,  as  shown  in  Fig.  17.  The  assumed  mean 

flow  deficit  profile  behind  the  struts  is  shown  in  Fig.  18.  The  wake 

w 

turbulence  convection  velocity  is  set  to  be  -  u^/2,  and  the 

turbulence  intensity  variation  according  to  the  distance  from  the 

v  2*  1/2 

trailing  edge  of  the  upstream  strut  is  given  by  ((u  )  )  /U^  * 

0.08(2x/e  +•  0.05)  *  ,  where  x  is  the  coordinate  along  the  airfoil 
measured  from  the  trailing  edge,  and  c  is  the  length  of  the  airfoil, 
see  Ref.  [44].  In  the  computation  shown  in  Figs.  14  through  16,  no 
consideration  is  made  of  the  wake  turbulence  length  scale  variation 
with  changes  in  rotor-strut  spacing  d.  In  fact,  the  lengch  scales, 

A^  and  A  ,  are  too  small  to  represent  a  blade  with  a  dipole;  indeed, 
the  dipole  domination  may  not  even  be  attained. 


Non-Dimensionalized  ^S^,(r,$,t  ,u))  for  B  -  3  and 


„/ A, ft  =  0.0542 


-Dependence  of  Non-Diinens 


Non-Dimensional  ized  <SfS>  (r  ,4  ,  t  ,w)>  for  B  =  9,  and 


where 


a/2 


A  -  0.68c{C  (x/c-0.35)}  '  /S8 

CQ  =  airfoil  profile  drag  coefficient 

c  =  airfoil  chord  length 

x  =  coordinate  along  airfoil  measured  from  midchord 
Fig.  17.  Wake  Position  Function  Profile 
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IS.  Mean  Flow  Deficit  Profile 


CHAPTER  VII 


CONCLUSIONS 


The  time  averaged  intensity  density  function  of  the  acoustic 
radiation  from  rotating  blades  was  derived  by  replacing  blades  with 
rotating  dipoles.  The  following  conclusions  can  be  drawn  from  the 
derivation : 

1.  The  rectilinear  model  of  a  rotating  blade  row  is  applicable  only 
for  the  frequency  range  such  that  u>>>Q;  see  Eq.  (2.25). 

2.  The  effect  of  the  variation  of  the  distance  between  the  observer 
and  rotating  dipoles  can  be  neglected  when  R-.'sino/a  <<  1;  see 
Eq.  (2.26)  through  Eq.  (2.28).  However,  care  must  be  taken  as 
exhibited  in  Fig.  10. 

3.  The  retarded  time  difference  of  the  lift  among  blades  and  its 
derivative  can  be  neglected  when  the  dipole  source  dominates  over 
other  sources,  the  speed  of  the  rotor  is  low,  and/or  the  aspect 


ratio  is  low;  see  Eq .  (5.5)  and  (5.6). 

4.  Two  non-dimensional  parameters,  and  R  /A^,  govern  the 

proraie  of  the  noise  spectrum.  As  L’j/'l!  decreases,  the  level  of 
the  blade-passing  frequency  ana  its  multiples  increases,  and  as 


R^/.V  increases  the  level  also  increases;  see  rigs.  5,  6,  and  7. 

5.  The  non-dimensionalized  time  average  acoustic  intensity  density 

IT:  2  3 

function  is  proportional  to  ( fu  )"/a  )(R_.'./a)  under  the  condi¬ 


tion  that  other  non-dimensional  parameters  are  fixed; 


see 
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Eq.  (5.22). 

6.  The  effect  of  the  wake  turbulence  can  be  reduced  by  reducing  the 

w 

wake  width,  increasing  decreasing  R  /A^,  decreasing  the 

v  2  i  ? 

number  of  struts,  and  reducing  the  ratio,  (u  )  /(u  )“;  see 
Eq.  (5.22). 

7.  The  number  of  struts  and  rotor  blades  should  be  chosen  so  that 
the  least  common  multiple  of  S  and  3  multiplied  by  ,1  is  in  the 
cutoff  range  of  the  Sears'  function  in  order  to  minimize  the 
effect  of  the  mean  flow  deficit;  see  Eq.  (5.22). 

8.  The  acoustic  intensity  at  the  off-angle  from  the  rotor  axis  is 
not  stationary  even  if  the  inflow  turbulence  is  frozen  and 
homogeneous;  see  Eq.  (5.11).  The  acoustic  intensity  from  the 
rotor  with  struts  is  not  stationary  even  at  the  rotor  axis;  see 
Eq.  (5.22). 

9.  Increasing  B,  and/or  increasing  2  under  the  condition  that  L’  / A  2 

3  i 

is  fixed  increases  the  bandw  .th  of  the  peak  of  the  blade-passing 
frequency  and  its  multiples;  see  Feg.  8. 

10.  To  find  the  effect  of  the  rotor  strut  spacing  on  (r ,  } ,  t  ,w))  , 

consideration  of  the  wake  properties  and  mean  flow  deficit 
profile  with  respect  to  the  spacing  is  necessary. 

11.  As  was  derived  in  Eq.  (5.22),  (r , ? , t ,u)}  is  expressed  by  the 

terms  due  to  the  inflow  turbulence  and  wake  turbulence,  and  the 
term  due  to  the  mean  flow  deficit.  Hence,  we  can  deal  with  the 

(r,  J.t.-j))  due  to  the  inflow  turbulence,  the  wake  turbulence, 
and  mean  flow  deficit  separately,  although  this  separation  is 
hypothetical. 
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APPENDIX  I 


Given  Eq,  (2.17)  as  follows: 
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where 
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D  =  r  -  R  sino  cos(2t  +  i  )  , 
mn  n  m 


(AI.2) 


and  integrating  Eq.  (AI.l)  with  respect  to  t,  one  obtains 
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and  the  relation,  dim (?(t)) }/d'?  *  (dm/dr) (di/df) ,  is  used.  Consider¬ 
ing  the  terms  in  the  parenthesis  of  Eq.  (AI.3),  one  obtains  the 
following  equation: 
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where  D  is  approximated  by  r,  due  to  the  assumption  that  R  <<  r  for 
mn  n 

all  n  *  1,2,... ,N  and  because  the  time  history  of  the  pressure  is  not 
sought.  Taking  the  x-Fourier  transform  of  Eq.  (AI.5)  and  taking  the 
ratio  of  the  1st  and  3rd  terms  yields 
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.  3D 

where  G(u)  *  "^Kg  (x)}  and  A (oj)  ="3,{1/(1+ - r-— ) ; .  Now,  by 

mn  3  dT 

utilizing  a  Taylor  expansion  on  A,  one  obtains 
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The  above  equation  shows  that  the  steady  state  component  has  a 
magnitude  "1".  Therefore,  because  3D  /a3x  <<  1  (since  a  low  speed 
rotor  is  considered  here),  G(ui)*A(ui)  =  G(w)  for  all  u>.  Hence,  the 
absoluce  value  of  the  ratio  of  the  lSt  and  3rd  terms  becomes  u>r/a. 

rcJ 

If  ur/a  >>  1,  the  3  term  is  negligible.  The  condition  wr/a  >>  1 
defines  the  acoustic  far  field. 


APPENDIX  II 


Lowson  [7]  derived  Che  following  Fourier  coefficients  of  the 

amplitude  modulation  for  3rt^  and  4C^  terms  in  Eq.  (2.25),  under  the 

condition  !>  =0, 

m 


cosfl(t  -  D  /a) 
mn 


R  a 

fl  +  n-  sind  siml(t  -  D  /a)) 
v  a  mn  * 


-  =  i 
3  p=- 


R  5 2 

-jpJ  (p - sino) 

_ . _ 2_  a 


R  .1 
n 


■sm<? 


jP(nt  +  it/ 2)  .  (Aii.i) 


Taking  the  phase  J  ,  of  the  amplitude  modulation  in  Eq .  (2.25)  into 
account  yields 


cos{P.(t  -  D  /a)  +  >  } 
mn  m 


R  .1 

00  -jpJ  (p— ^-sino) 

—  y  3 


r  ^  r 

(l  +  — sind  simr.(c  -  D  /a)  +  7  })3  p=_“  -2-s 

'■a  mn  m  1  a 


in® 


jp(Qt  +  t/2  +  J>  )  • 
e  m 


(AH.  2) 


Now,  consider  the  Fourier  coefficient  of  the  other  amplitude 

modulation  in  Eq.  (2.25),  for  instance,  the  Fourier  coefficient  of 

(l  +  (R  .'./a)sin®  sin..l(t  -  D  /a)  >  :)  “.  Taking  the  time 

n  an  m  '  ’ 

derivative  of  I  1  (R  3/a)sin®  sin-.P.(t  -  D  /a)  +  }  ;1  3  vields 

n  an  m  1 


.  ,  D'R  sin®  siny 

57' 1  +  (^a-Va)sino  sinv}"  - - ^73 -  . 

al!  +-^inO  siny)3  (An>3) 


where  the  following  relations  are  used: 
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d _ dy£T£_ 

dt  “  dT  dt  dy  ’ 


(All. 4) 


Y  -  n(t  -  D  /a)  +  i  . 

mn  m  ’ 


and 


T  =>  S2(t  -  r/a)  +  p  = 
m 


R  SI 

Y - sinpcosy 


(All. 5) 


(All. 6) 


By  comparing  Eq.  (All. 2)  with  Eq.  (All. 3),  one  obtains  the 
following  expansion: 


R  Q 

(i  +  -? 


ps»— oo  **  “ 


R  2 


Jp(Qt  +  tt/2  +  J  ) 


(All. 7) 


where  the  steady  state  component  can  be  computed  by 
o  2tt/Q  r  n 

JL  f  /i  ,  n  .  _ 7 

2t  J  (1  +  “7"  Sln®  siny)  Adt  =  1  , 

R  SI 

since  dy  «  (1  +  -~-sind  siny)-1dt. 

Therefore,  one  obtains 


(All. 8) 
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(All. 9) 
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Now,  consider  the  Fourier  coefficient  of  sin(i2(t-D  /a)+‘t>  }/ 
RnQ  2  nn  m 

(l  +  — ^-sin?  s±n{2  (t  -  D^/a)  +  $  })  •  Since  it  can  be  expressed  by 


sin{0(t-D  /a)-i-5  } 
mn  m 


R  3  R  ?i 

1  +  sinosin{i2(t-D  / a)-K>  }  1  +  — sinpsin{£l(t-D  /a)4-t>  } 
a  mn  m  a  mn  m 


the  Fourier  coefficient  can  be  obtained  through  coefficients  for  each 
term  in  the  above. 

R  Q 

Taking  the  time  derivative  of  sin{0(t-D  /a)+i>  } /  f  1  +  n  sinp 

mn  m  11  a 


sin{H(t-D  /a)-K>  })  yields 
mn  m  1  J 


smv 


^cosv 


dt 1  R 


n  RnP‘  3 

1  +  — 2— sinosiny  (l  +  — sindsiny) 


(All. 10) 


where  Eq.  (All. 4),  Eq.  (All. 5),  and  Eq.  (All. 6)  are  utilized. 

By  comparing  Eq.  (All. 10)  with  Eq.  (All. 2),  the  following 
expansion  can  be  obtained: 


smy 
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1  +  - smbsiny  r  - smO 
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(50p  -  1)  ,  (All. 11) 


where  the  steady  state  component  is  obtained  by 


2-r/Q 


smv 
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(All. 12) 


Therefore,  utilizing  Eq.  (All. 9)  and  Eq .  (All. 12),  the 


following  expansion  is  obtained: 
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APPENDIX  III 


The  expected  value  of  Rect(--*  -)  is  computed  by  the  following 
integration: 

oo 

g(t)  =  //  Rect k)p  (b)p  (w)dbdw  , 

_«>  w  r  r 

where 

Pr(b>  -  Ibi  <  A 

.0  otherwise  , 

L  w  L  a2 

otherwise  ,  (AIII.3) 

and  b  and  w  are  random  variables  for  the  position  and  width  of 
Rect(t),  respectively. 

Integrating  with  respect  to  b  and  w,  one  obtains 


(AIII.l) 


(AIII.2) 
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Now,  consider  the  special  case  such  that  p^Cb)  =  6(b)  and 


6(w  -  «c) .  Then,  g(t)  becomes 


g(t)  *  Rect(t/wc)  . 


(AIII.5) 


Eq.  (AIII.5)  represents  the  case  where  Rect(t)  is  a  deterministic 


function  whose  width  and  position  are  2w^  and  zero  in  the  time  axis, 


respectively. 
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